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We give the optimal bounds on the phase-estimation precision for mixed Gaussian states in 
the single-copy and many-copy regimes. Specifically, we focus on displaced thermal and squeezed 
thermal states. We find that while for displaced thermal states an increase in temperature reduces 
the estimation fidelity, for squeezed thermal states a larger temperature can enhance the estimation 
fidelity. The many-copy optimal bounds are compared with the minimum variance achieved by 
three important single-shot measurement strategies. We show that the single-copy canonical phase 
measurement does not always attain the optimal bounds in the many-copy scenario. Adaptive 
homodyning schemes do attain the bounds for displaced thermal states, but for squeezed states 
they yield fidelities that are insensitive to temperature variations and are, therefore, sub-optimal. 
Finally, we find that heterodyne measurements perform very poorly for pure states but can attain 
the optimal bounds for sufficiently mixed states. We apply our results to investigate the influence of 
losses in an optical metrology experiment. In the presence of losses squeezed states cease to provide 
Heisenberg limited precision and their performance is close to that of coherent states with the same 
mean photon number. 
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I. INTRODUCTION 

Since the early work of Dirac [l[ the problem of mea- 
suring the phase imprinted on a state of light has been 
a subject of great debate. The problem was to find the 
quantum counterpart to the classical phase observable. 
The observable should be the conjugate variable of the 
number operator n, and the corresponding uncertainty 
relation could be derived. From a quantum measure- 
ment point of view the lack of such a self-adjoint observ- 
able was overcome with a different approach based on 
quantum estimation theory [H, El, GJ , which gives a pre- 
cise and operational meaning to "measuring" the phase 
(without resorting to the notion of observable) and pro- 
vides reasonable uncertainty relations. On top of this 
fundamental motivation, phase estimation is at the heart 
of many quantum metrology Q applications, such as im- 
provement of frequency standards [6| , gravitational- wave 
detectors 0,11], and clock synchronization |§ [Toj . and is 
strongly related to quantum computation [llj and quan- 
tum cryptography 12| . It is therefore an essential task to 
compute the efficiency and ultimate bounds on the pre- 
cision of phase-estimation. Optical implementations of 
quantum metrology applications are within reach of cur- 
rent technolo gy (see, for ex amp le, recent experimental 
achievements pi, QJ, QJ, GJ, df ) • 

Here we study optimal protocols to estimate the phase 
encoded in states of light. In particular, we will deal 
with two families of continuous variable (CV) states that 
arc described by a Gaussian characteristic function: dis- 
placed thermal states and squeezed thermal sates. Gaus- 
sian states are mathematically easy to handle and very 
relevant experimentally because they provide a very good 
description of the states available in the laboratories: 
laser fields manipulated with passive and active linear 
optical elements. Although there is some previous work 



in the case of pure states (see [18| and references therein) , 
very little is known about estimation in mixed states, 
which are present in realistic unavoidably noisy scenar- 
ios. This is the issue we tackle in this work. 

In the next section we present the general framework, 
introduce our figures of merit, and show that the max- 
imum value of the average fidelity can be achieved by 
the single-seed covariant generalized measurement. In 
Sec. IIIII we apply the techniques introduced in Sec. [TT] 
to estimate the phase encoded in a single copy of a dis- 
placed thermal state (also known as coherent thermal 
states). In Sec. II VI we proceed by analyzing the problem 
for squeezed thermal states of a specific class: those that 
emerge when a squeezed vacuum state is sent through 
a lossy channel. Contrary to intuition and in contrast 
with the case of coherent thermal states, we find that 
for squeezed thermal states the estimation precision im- 
proves with the temperature of the encoding state. In 
Sec. [V] we compute the maximum fidelity when many 
copies of the state are available and study the perfor- 
mance of three particular single-shot strategies that only 
use individual measurements on each copy. Finally in 
Sec. I VII we use our results to evaluate the effect of losses 
in metrology experiments. We end with the main con- 
clusions of our work. 



II. GENERAL FRAMEWORK 

In this section we introduce the general problem of 
estimating a phase. In our setting a system evolves under 
a unitary transformation 



P {<j>) = u{4>) P u{4>) 



(i) 



where p is a general Gaussian state and U (4>) is the uni- 
tary operator U((j>) = e^ n , with <f> € [0, 2ir) and n = a^a 



2 



is the number operator, and a) and a are the creation and 
annihilation operators that satisfy the bosonic commuta- 
tion relations [a, at] = 1 and [a, a] = [a^a^] = 0. Al- 
though our focus here is on phase-shifted Gaussian states, 
the results in this section can be applied also to non- 
Gaussian states, in particular also to finite dimensional 
states, i.e. generic qudits [see however the comment after 
equation JSJ]. 

Information on the phase is obtained through a mea- 
surement specified by a positive operator- valued measure 
(POVM), i.e. by a set of operators {O x } that are pos- 
itive semi-definite and add up to the identity. The out- 
come probabilities depend on the measured state through 
Born's rule p{x\p) = tr(O x p). Estimation is thus inher- 
ently not perfect, and we need a figure of merit to quan- 
tify how close is our guess (f) x , based on outcome x-i to 
the real value of <j>. We will take the usual function 

fl{(j),<j) x ) = cos[l((j}-4> x )), (2) 

where I = 1 for displaced states and I = 2 for squeezed 
states. The factor of 2 in the second case takes into 
account the symmetry of squeezed states under a phase- 
shift of 7r, due to the their double-photon structure (see 
below). We will loosely refer to /; as the estimation fi- 
delity. The corresponding average fidelity is given by 

*i=£ f n ^m^tTwm (3) 

Next we proof that the maximum average fidelity can 
be attained with a covariant POVM and reads as 

oo 

^ = X>«.«+'I- ( 4 ) 
n=0 

In fact we can derive the maximum average fidelity for 
a rather general family of merit functions: f{4 > i4>x) = 
J2i aifi{4>, <t> x ) = J2i a i cos[l((f> - <f> x )], with a t > 0. This, 
of course, includes the above choices as well as most of 
the commonly used figures of merit. To this end we write 
the input state and POVM elements in the Fock basis, 

p(4>) = J2pn,n^ (n ~ n ' ) \n}(n'\, (5) 

O*, n = (n'\O x \n), and develop the expression for the 
average fidelity to arrive to the following upper bound 

X-,n,n' 

= E /^i«»[I(0-^)]e^-"'W^,„ 

X-,n,n' J 
X-,n,n' ,1 

<Y, ai \p n , n+l \\(% +l J. (6) 

l,n,x 



The inequality is saturated if and only if the following 
relation between the phases 0^ n = arg n and VVn,n = 
arg p m<n holds: 

^n,m = V'n.m + ^ x (n. - m) Vn, m. (7) 

Now, from the positivity of O x and the fact that the ge- 
ometric mean is bounded from above by the arithmetic 

mean, it follows |0*J < y/0* n O*, m < 1/2(0* n + 
Om m )i which together with the POVM completeness re- 
lation J2 X On, m = <^n,m leads to 

^^E ai Ei^"+'i- ( 8 ) 

l n 

This bound can be attained iff 1 0*^ „| = 1 Q* m \ = \ O x n \ 
for all to, n, and x- This is precisely satisfied by the well- 
known (canonical) phase-measurement Q. This is given 
by a continuous POVM [labeled by 9 <E [0, 2ir) instead 
of a discrete index x] with elements (f = l/(2ir) \6)(9\ 
where \9) = J2 n e%dn \ n )- Condition (J7]) on the phases 
is automatically met if the density matrix p has posi- 
tive entries in the Fock basis, i.e., ip n ,m = 0, and the 
phase 9 is guessed whenever the measurement outcome 
Cf occurs. If p does have non-zero phases of the form 
i>n,m — C( n ) ~ C( to )j i- e - they can be removed by a 
unitary transformation U = ^2 n e 1 ^") \n)(n\, the phase 
measurement given above needs to be corrected by the 
corresponding unitary, i.e., \8) = J2 n e t '- ne+ ^ n '' \n). For 
more general phases of p n . m , the bound holds if condition 
(J7J) does not conflict with the positivity condition on the 
POVM elements. For any dimension d > 2 one can easily 
find examples of seed states p(0) with arbitrary phases 
for which the bound cannot be attained. However, for the 
Gaussian states under study here, phases can always be 
unitarily cancelled, and therefore ^ does indeed give the 
maximum fidelity. In fact, for the figures of merit that 
we use in this paper, and in general for those that only 
have a single Fourier component (only one non-zero co- 
efficient ai), the bound can always be attained since any 
arbitrary phase ip„_ n+ i can be unitarily reabsorbed by the 
POVM, i.e., \0) = J2 n e^ n0+ ^ \n), where the phases 
satisfy the recursion relation £(n + I) = £(n) + ip n , n +i, 
with £(n) = for n < I. 

Note that the phase measurement is optimal not only 
for a large family of input states but also for a very 
wide class of figures of merit. A remarkable excep- 
tion occurs for the state estimation fidelity f(p x ,P4>) — 
{^A^/P~x\ r P<i>\) 2 > where p x is the state "guessed" on mea- 
surement outcome x- This optimal guess state is not nec- 
essarily one of the possible input states and the figure 
of merit cannot always be written as a Fourier series with 
positive coefficients. Therefore the above derivation does 
not hold. In fact, in state estimation of phasc-covariant 
states there are examples where the phase measurement 
is known to be sub-optimal [l9l |. 

Equation @ provides the ultimate bounds in phase 
estimation; however the phase measurement saturating 
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this bound may be very difficult (if not impossible) to im- 
plement experimentally in a single-shot (or single copy) 
measurement [2(|. However, as we will sec in section f7\ 
when a large number of copies are available there ex- 
ist known measurement schemes that attain the optimal 
bounds in some particular cases. 

The remaining of this section is devoted to Gaussian 
states of light. The Gaussian states are those that are 
fully characterized by the first and second moments of 
the field quadratures, i.e. by the displacement vector d = 
tr(Rp) and the covariance matrix Ty = tr(p{i?j.— df., Ri — 
where R 1 = l/y/2{a + a)) and R 2 = i/V2(al - a) 
are (conjugated) quadratures. 

An equivalent, perhaps more operational, definition 
can be given in terms of the action of the squeezing op- 
erator S(r) = exp[|(a 2 — a* )], and the displacement 



operator D(a) 



~ a a on a thermal state 



(9) 



namely, 

PftvM = U(<f>)D(a)S(r)pf)S(ry D(a^U(^ . (10) 

For a general Gaussian state the mean photon number 
can be easily seen to be 



(n) = H 2 + np + {2n l3 + 1) sinh 2 r, 



(11) 



where np = (e* 3 — 1) 1 is the thermal mean photon num- 
ber. The covariance matrix of a thermal state is simply 
Tp =701, with 70 = tanh~ 1 (/3/2) = 2710 + 1. The 
squeezing operator S(r) decreases the fluctuations in one 
quadrature and increases them by the same factor in the 
second quadrature: 



r ri/3 =(2710 + 1) 







,2c 



(12) 



The displacement D(a) does not change the second 
moments, instead it only induces a displacement d = 
V / 2(Rca, Ima) T . Finally the phase operator U(<p) pro- 
duces a rotation in phase space which induces the corre- 
sponding transformation on the displacement vector and 
covariance matrix : d' = O^d and T' = O^TOT with 



O, 



COS ( 

sin t 



— sm< 
cos 6 



(13) 



We will next concentrate on two classes of Gaussian 
mixed states: displaced thermal states, also called coher- 
ent thermal states, and squeezed thermal states. 



III. COHERENT THERMAL STATES 

Coherent states, defined as 

2 OO n 

|a)=I)( a )|0>=e-^£4=|n>, (14) 



are very relevant in CV implementations of quantum in- 
formation protocols since they provide a very good de- 
scription of the states produced by a laser: they are 
states with a well-defined amplitude and phase, and 
with minimal fluctuations in both quadratures (see how- 
ever [2l| and [lij for some caveats on this description 
of the output of a laser). A displaced thermal state 
Pp,a = D{a)ppD{a)^ is also well represented by an am- 
plitude and a phase but fluctuations are larger. 

In order to calculate the maximum fidelity ^ it is 
convenient to write the thermal state in the P-function 
representation 



1 

imp 



= — I d 2 a'e ~ \a'){a'\ 



(15) 



One can check that this is indeed a thermal state just by 
computing the matrix elements in the Fock basis 



1 f i°'i 2 
(n\pp\m) = / d ae (n\a)(a\m) 

7T710 J 



(16) 



" w "(l + n /J )"+i" 
In the Fock basis the coherent thermal state reads as 

1 f i°'i 2 
PBa = d 2 a'e n P W + a)(af + a\ = (17) 

7T710 J 



E 

k.i 



d 2 a' c -^ c - W+a] * (a' + a) V + a)** 

7T710 \fk\\/T\ 



Equation ([J), for 1 = 1, gives the maximum fidelity, which 
reads: 



imp 



fc!VF+T 



(18) 

At this point we use the following integral representation, 
1 1 



Vk + 1 



Vt 



(19) 



that allows us to evaluate the sum. The integral over 
a' can then be easily evaluated: the imaginary part is 
trivially zero and the real part is a Gaussian integral 
that is equal to 



dy 



^/-Hl-y/(l-npy)} 



, (20) 



after the change in variable y = (1— e _ *)/[l+7ig(l— e - *)]. 
This expression can be evaluated numerically with ar- 
bitrary precision for different thermal photon numbers 
np. Figure [1] shows the average fidelity as a function of 
the field amplitude for different temperatures. As one 
could expect, higher temperatures increase the fluctua- 
tions without changing the field amplitudes and therefore 
decrease the precision of the estimation. 
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FIG. 1: (Color online) Average fidelity for coherent thermal 
states for different np: np = 0.1 (solid), np = 0.5 (dashed), 
np — 1.5 (dotted), np — 2.5 (dash-dotted). 



The asymptotic behavior of the fidelity at large field 
amplitudes can be computed analytically from (|20|) by 
expanding the non-exponential factor of the integrand 
up to second order, 



T x ~ 1 



2np + 1 
8n a 



(21) 



where n a = \a\ 2 is the mean photon number of a coherent 
state | a), or the contribution of the displacement D(a) 
to the mean photon number of a coherent thermal state 
n = rip + \a\ . In the limit of zero temperature, np = 0, 
we recover the result for pure states T\ ~ 1 — l/(8n) 
derived in 23]. 

We next turn to the opposite regime. In order to com- 
pute the fidelity for low- field amplitudes (a — > 0), we 
write the thermal state in the Fock basis and expand 
D(a) up to first order in a, D(a) w 1 — a(a — af). With 
this one can easily obtain a closed expression for the fi- 
delity up to first order in \a\. For simplicity here we only 
give its asymptotic behavior at high and low tempera- 
tures: 



7T / n n 



2 V 2np + 1 



for np > 1 



(22) 



i a [l - (2 - V2)n p ] for np < 1. 



IV. SQUEEZED THERMAL STATES 



Squeezed vacuum states are defined as 

oo 

|r)=5(r)|0) = (l-A 2 ) 1 / 4 ^ 



n=0 



\2n), (23) 



where A is related to the squeezing parameter r and the 
mean photon number n r = sinh 2 r by A = tanh r = 
n r j (n r + 1). Squeezed states are those for which the 
uncertainty (or fluctuations) on a given quadrature is re- 
duced below the standard quantum limit (the quantum 



fluctuations of vacuum) at the expense of increasing the 
uncertainty in the conjugate quadrature. One can there- 
fore expect an enhanced performance of squeezed states 
in phase estimation or in other high-precision applica- 
tions [1, [24|, HH . This states exhibit very non-classical 
features, such as sub-Poissonian statistics or the ability 
to generate entanglement [26|, [27| , which allows for a va- 
riety of applications in quantum information theory [28| . 
Squeezed states have been successf ully produced in the 
laboratories for some time now [29|, [3(| reaching squeez- 
ing of up to 10 dB (corresponding to r < 1.15 or n r < 2) 



171 ] . Squeezed thermal 



in current experiments [151 . \lt 
states are those obtained by squeezing an initial thermal 
state (instead of the vacuum), 



pp r = S(r)ppS(r)K 



(24) 



Therefore one expects quantum features to be less pro- 
nounced. At high enough temperatures [for e~ 2r (2np + 
1) > 1] the state will even become classical [3l|, i.e. a 
mixture of coherent states, loosing its entanglement po- 
tential [321 ] and other quantum features. 

In order to obtain the maximum fidelity ([5]) for 
squeezed thermal state we only need to sum the outer 
diagonal elements |p„,„+2|- A closed form expression for 
the matrix elements (m\ S(r) \n) can be found in (33j for 
arbitrary Fock states, but it is highly non-trivial and 
turns the evaluation of the fidelity into a very hard com- 
putation. Let us hence consider a much more tractable, 
but still very relevant, family of squeezed thermal states. 

Most sources of noise in quantum optical experiments 
(e.g., mode mismatch, misalignment, absorption of opti- 
cal elements, and non-unit detector efficiencies) can be 
cast in terms of linear losses. The family that we will 
study here consists of the mixed states that arise when 
a squeezed vacuum undergoes losses during the course of 
an experiment. Since linear losses can be modeled by a 
beam splitter (BS) of transmittance T, this family has 
the following simple characterization: 

where B = exp[9(a^b - a6 f )], with T = 1 — R = cos 2 6, 
is the BS transformation that acts on both the system 
(mode a) and the environment (mode b, with [b, w\ = 1), 
which is initially in the vacuum state. Notice that the 
effect of losses commutes with the phase operation and, 
therefore, the states under consideration are of form (fTJ), 
i.e., p((p) = U (<p)pU ((/))' . In what follows we find how the 
input squeezing parameter ro, and the transmittance of 
the channel, T, are related to the inverse temperature, /3, 
and squeezing parameter r of the output state. 

The covariance matrix associated with the two-mode 
initial state, |r ) a ® |0) h , is given by r in = r rOi0 © 1. The 
action of the BS transforms the quadratures according to 
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the symplectic operation, 



leading to the two-mode covariance matrix 



Vo 



cos 9 





sin 9 








cos 9 





sin 


— sm9 





cos 9 








— sin 9 





cos 



(25) 



v e r in v e 



/ e- 2ro T + R (e~ 2r ° - 1)VtR _\ 

e 2r "T + R {e 2r °-l)VTR 



(e-' 2r « - 1)VTR 

\ o 







{e 2ro - 1)VTR 



e- 2r( >R + T 






ro R + T J 



(26) 



If we consider only mode a we obtain 



1 out — 



e -2r oT . 





R 

e 2r °T- 



i? 



(27) 



That is, the effect of losses turns the initial squeezed vac- 
uum |ro) into a squeezed thermal state pp, r . Equating the 
covariance matrix (|27|) to its general form for squeezed 
thermal states (fl2"]) we find the following equations relat- 
ing input and output parameters, 



2ng + 1 = tanh 



-10 



y/T 2 + 2TRcosh2r + R? , 



,2r 



e 2r «T + R 
e- 2r °T + R' 



(28) 



Note that equations (f2"8")) do not have solution for all val- 
ues of np and r. That is, not all squeezed thermal states 
can be viewed as a squeezed vacuum that has degraded 
through a lossy channel. 

Now we can write the density matrix in the Fock basis 
for this family of squeezed thermal states by sending a 
squeezed vacuum \r ) ([125]) through a BS (f2"5|) and tracing 
out mode b, 



PpA<t>) = V 1 ~ 



*8 E (-T 



n-j-m 



(29) 



ilm.i 



J2\/PT n (k)JpT m (k) \2n-k}(2m-k\ 



with Ao = tanh ro and with the binomial distributions 
J5y"(fc) and p 2 ^ n (k) given by 



p£(fc) = 



AH 



(2V-fc)!fc! 

Using the maximum fidelity can be written as: 

2n+l 



(30) 



1 



A oE 



(2n)!(2n + 2)! R k r( 2n+1 " fe ) 



n!(n + l)!fc!(2n - fc)V2n + 2 - ky/2n + 1 - fc 



(31) 



At this point, we use the integral representation (|19p for 
both l/y/2n + 2 — fc and l/y2n + 1 — fc, which reads as 

1 1 



\/2n + 2 - fc V2n + 1 - fc 

I f°° ^fi e ~x(2n+2-k) f°° e _j,(2n+l-fe) > ^32) 
K Jq V% Jo VV 

This enables us to perform the sum over n and fc explic- 
itly. After the change of variables x — ut and y = (1 — u)t 
one easily identifies the integral representation of the 
Bessel function 



Io(t) = 



1 



du 



t(l-2u) 



7T Jo V^(l _ M ) 

and the fidelity can be cast in the compact form 

dte-^I (t/2) T 
{1-A 2 [i? + Te-*] 2 } 3 / 2 ' 



•7"2 = Aq y 1 — A 



(33) 



From Equations (|28p relating r and with tq and T it 
is straightforward to calculate the fidelity for states that 
have the same squeezing parameter r but different tem- 
perature. In Figurc[5]we show the fidelity as a function of 
r for a pure state, |r), and for a mixed state. The latter is 
taken to be the state that results from sending |ro = 2?') 
through a 50:50 BS (T = 1/2), which indeed has a squeez- 
ing parameter equal to r. We see that the mixed states 
(non-zero temperature) allow for better phase estimation 
than pure states (zero temperature). This behavior is a 
bit puzzling at first sight and opposite to what we found 
for coherent states. It is true that noise per se increases 
the uncertainty in any estimated quantity. However, in- 
creasing the temperature has a second effect that goes 
in the opposite direction. When the squeezing operator 
acts on a thermal state the number of photons increases 
in a non-linear fashion, 



(n) = rip + (2np + 1) sinh 2 r. 



(34) 



This boosted increment gives rise to an improvement on 
the phase sensitivity that outweighs the adverse effects 
of the temperature. In the case of coherent states, the 
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action of displacement adds the same amount of energy 
independently of the temperature, and hence the detri- 
mental effect of thermal noise is not counterbalanced 
by the energy increase. Figure [5] also shows (dotted 
line) the fidelity of a pure squeezed |f) state with the 
same mean energy of the squeezed thermal state, i.e., 
rif. = sinh 2 f = 1/2 (2r| n \2r) = 1/2 sinh 2 (2r), where we 
have used that the BS is balanced and half of the pho- 
tons are lost. We find that, given a fixed mean energy, 
it is clearly more advantageous to increase the squeezing 
parameter rather than the temperature. 




FIG. 2: (Color online) Fidelity for pure squeezed state \r) 
(solid) and for states with the same squeezing parameter 
emerging from a 50/50 BS (dashed). We also show the fi- 
delity for a pure squeezed state \f) with the same energy as 
the former thermal states, i.e. (n) = sinh 2 r = 1/2 sinh 2 2r 
(dotted line). 

The fidelity behavior in the asymptotic limit when the 
squeezing parameter is large can be computed analyti- 
cally. With the change in variables t = — law, Eq. 
becomes 



J~2 — -Vn/l 



a 2 t 



dw ^/w \ m w ) 



Ag(iZ- 



TV) 2 ] 3 / 2 ' 



(35) 



Note that in the limit we are interested in, Ao approaches 
unity and the dominant contribution to the fidelity comes 
from the region u> ~ 1. It is, therefore, convenient to sep- 
arate this contribution by writing the numerator of (|35|) 
as 1 — [1 — \/wIq(^ law)] and integrating the two terms 
separately. The integral corresponding to the first term 
(unity) is straightforward and contributes to both lead- 
ing and subleading orders. The term in square brackets 
goes as 1 — w for id w 1 and is subleading; thus we can 
safely put Ao = 1 in the corresponding integral, as the 
neglected terms will give even higher order contributions 
to JT 2 . The result can be cast as 



^2 = 1- 



where 



R 



-T 



1 - s/wl {\ \tlw) 



dw , 

o [i - (R 



•TVj 2 ] 3 / 2 



(36) 



(37) 



The remaining integral can be computed numerically 
to arbitrary accuracy for any value of the channel trans- 
mittance T. For T = 1 (pure squeezed state) we ob- 
tain £(1) = 0.55, in agreement with [23|. We further 
note that, as expected, £ is a decreasing function of T. 
Thus, e.g., £(1/2) = 1.20 and £(1/3) = 1.62. By inter- 
polating £(T) between the two regimes where it can be 
analytically computed (T — ► 1 and T — > 0) one can check 
that to an accuracy of more than 98%, one has 



R 



2T 



C2T, 



(38) 




where C\ = 0.54 and c 2 = 0.17. 

Let us now consider the opposite limit of low energy, 
which is relevant because the amount of squeezing avail- 
able in laboratories is often quite limited. For Ao <C 1 we 
can take Ao ~ \/no and retain only the first term (n = 0) 
in the sum (|31~j) to get 



(39) 



or T<x w Tro/v2 f=a r/y/2 in terms of the squeezing 
parameter. This explains the linear behavior shown in 
Fig- HI f° r small values of r. 



V. MANY COPIES 

The calculation of the fidelity when several copies N 
are available is in general a hard task. The case of pure 
states has only been solved recently [23| by computing 
the matrix elements of in a nondegenerate eigen- 

basis of n t — y^-i <x|di- The case of mixed states be- 
comes further involved since the basis that spans p® N 
is degenerate with respect to n t . Here we will follow a 
different approach to calculate the asymptotic behavior 
of the fidelity in the limit of large number of copies. 

We first notice that in this limit the estimation be- 
comes very accurate (</> — <p x will be typically small) and 
therefore we can relate the estimation fidelity ([!]) to the 
variance: T\ = 1 — Var[0]/2 for coherent states, and 
T2 = 1 — 2Var[</>] for squeezed states. The Cramer- 
Rao bound [U, a well-known result in classical statis- 
tics, gives a lower bound on the variance of an unbiased 
estimator in terms of the Fisher information /(</>), which 
is a functional of the parametric family of probability 
distributions pixl^) from where the samples are drawn 

m=jd XV{ m(^^-)\ (40) 



For a large number of samples N and under some regu- 
larity conditions, the Cramer-Rao bound is attained and 
we have that Var opt [<j)\ = (NI) -1 . In their seminal pa- 
per Braunstein and Caves [35[ made use of this result 
to proof that if N copies of a state p(4>) are available, 
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the optimal estimation of <p based on the most general 
quantum measurement is 



Var°P t [0] = d^/(4iVd4 tJ ), 



(41) 



where the Bures metric ds^y can be obtained from 
the Bures distance between two infinitesimally close 
states: ds BU = 1 — F(p,p — dp), where F(pi 7 p 2 ) = 

tr {\ZVPiP 2 VPz) 1S tne q uan tum fidelity [36|. Closed 
expressions for this fidelity have been computed in 
[37L |H, [39| for general Gaussian states. Therefore, in 
order to obtain the maximum estimation fidelity we only 
need to compute F[p(0), p(d<j>)] for infinitesimally close 
Gaussian states. 

For displaced thermal states, pp, a , we find 



1 



1 



1 2np + l 



AN |a| 2 tanh/3/2 AN n a 



(42) 



For pure coherent states we recover the known results 
[l8| which agree with the standard quantum limit or 
shot noise limit (Var[</>] ~ We notice here that 

for pure states the iV-copy result follows straightfor- 
wardly from the single-copy case since an TV-copy co- 
herent state \a)® N is unitarily equivalent to a large 
amplitude coherent state \y/Na). The unitary can be 
simply implemented by a linear multiport device that 
realizes a linear mode transformation such that ao = 
l/y/~N(ai + <22 + . . . ajv)- For mixed states the error in- 
creases linearly with the size of the quadrature fluctua- 
tions, i.e. the phase variance scales linearly with the num- 
ber of thermal photons np. We also notice that although 
the above strict equivalence between the single-copy and 
many-copy regimes does not hold for mixed states, the 
variance at large field amplitudes has exactly the same 
behavior in both regimes [compare (|4"2"| with (|2ip]. 
For squeezed thermal states, pp, r , we find, 



Var°f 



[<t>] 



1 



1 



2N (1 + cosrT 1 /3) sinh 2 2r 



1 



l6Nn r (n r + 1) 



1 



{2np + l) 2 



(43) 



where we used that sinh 2 (2r) = An r (n r + 1). A gain , for 
pure states we recover the well-known result in [18], [If]] , 
which has the scaling of the so-called Heisenberg limit 



Var[0] 



Here, as in the single-copy scenario, we 



find that even though the phase-space distribution clearly 
shows that temperature increases the size of the fluctua- 
tions, a higher temperature prior to squeezing does in fact 
improve the estimation fidelity Notice that this behavior 
appears at all temperatures and persists even for classical 
states, i.e., when (2np + 1) > exp(2r) (see above). 

Finally, we would like to write the minimal variance in 
terms of the parameters of the lossy channel described 
in the previous section, i.e., the squeezing parameter r$ 
of the initial squeezed vacuum and the transmittance T 
of the channel. Using relation (|2"8"j) the optimal variance 



can be written after some algebra as 

Va^M= 1 + 2r(1 -^ Sinh2r ° ! (44) 
sq m 2 /vT 2 sinh 2 (2r ) ' V ' 

or in terms of the input's mean photon number hq = 
sinh 2 tq , 



VaxM 



1 + 2T(1 - T)n 



8NT 2 n (n + l)' 
which in the limits of high and low squeezing gives 



(45) 



1 



VaxM = 



8NT 2 n 



for n > 1, (46) 
for n < 1. (47) 



This shows that the Heisenberg limited precision cannot 
be attained in the presence of losses. 

The results we have presented so far for the many- 
copy scenario are again theoretical bounds on the optimal 
variance based on relation (|4"Tj) . The proof of this relation 
is constructive in the sense that it provides a particular 
single-copy measurement and an estimator that saturate 
the bound. However, this optimal measurement is given 
in terms of an observable called symmetric logarithmic 
derivative [35j . which typically cannot be implemented 
experimentally. In the following subsections we study 
three particular single-shot measurements and compare 
their performance with the optimal bounds. 



A. Heterodyne measurements 

Heterodyning is one of the first and most common de- 
tection schemes for state reconstruction in CV systems 
for it yields direct information about the phase-space dis- 
tribution (Q-function) that is enough to completely re- 
construct the state. The POVM of an ideal heterodyne 
measurement is given by { 1/ir | a)(a | } . 

The overlap between two arbitrary Gaussian states pA 
and pb is given by, 



tv(p APB ) = 2[dct{T A + T B )Y 



-s'crA+rs)- 1 * 



(48) 



where S is the difference between the displacement vec- 
tors of the two states 6 = dA~d,B, and T a and Tb are the 
covariance matrices of each state. Therefore, a hetero- 
dyne measurement on a squeezed thermal state pp,r{4>) 
will yield outcome a! . or more precisely its real and imag- 
inary parts {a' xl a' p } = {Rc[a'], Im[a']}, with probability, 

1 -d'evr-^d 
p(a'|# = -tr[|a>'| Pfi>r ] = j=^, (49) 

where d = \/2{a' x , a' p }, is defined in Q13p . and 



T, = 



V- 

v+ 



(50) 
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withV± = l + (2n + l)e ±2r . The Fisher informat ion can 
be computed from its definition in Eq. (|4T))) by taking the 
logarithmic derivative of the Gaussian (|49|) . 



da 



'2 e 6 [fl.utn.r-in J\l2 



[d^O^O^d)] 2 . (51) 



The derivative in the integrand is given by 

— =-iO <t> [o y ,T s jO = I — -— \ O+oJJ^ 

where we have used = exp(— i4>a y ) and a x and a y 
are the Pauli operators. Making the change in variable 
a' — ► a'e 1 ^, which corresponds to d — > O^e?, equation 
([S"Tj) becomes 



.„ 16a' 2 a' 2 
da' 2 21 p 



1 

V+) 2 4n r (n r 



1 

f l)(2n/3 



n r (2n/3 + 1) + (1 +n/3)" 



(52) 



For pure states the phase variance is given by 
1 



ANn r 



for np = (53) 



which has a shot-noise 0(n _1 ) scaling, instead of the 
Heisenberg 0(n~ 2 ) scaling that a squeezed vacuum state 
can provide: (|4"3")l . However, from (fB"2")) we also notice 
that, as in the case of the optimal measurement, the 
Fisher information increases with temperature. In par- 
ticular, in the limit of very high temperature we find, 



Var^M 



l6Nn r (n r + 1) 



for np — > oo, (54) 



which coincides with the optimal bound (|43p . So, het- 
erodyning performs very poorly (sub-optimal scaling) for 
high-purity states; but quite surprisingly it attains the 
optimal bounds at high temperature. 

Using (f2"5|) and (|52")) we can give the phase variance for 
an initial pure squeezed state of mean photon number no 
that is sent through a lossy channel of transmittance T: 



Va< et [0] = 



1 



(l-R 2 )n + 
ANT 2 n (n Q + 1)' 



(55) 



where we recall that R = 1 — T. 

In the limits of high and low initial squeezing we find, 



Var s h q et [</>] = 



2-T 



for no > 1, (56) 



ANTn 

Var sq 6t M = A Jr2 fOT n o « !• (57) 

Comparing with the optimal results in these regimes, (|46p 
and (|47|) , we see that in the low squeezing regime hetero- 
dyning performs a factor 2 worse than the optimal strat- 
egy, but the difference becomes even more pronounced 
for high squeezing. 



To calculate the heterodyne variance for coherent 
states we proceed along the same lines. The outcome 
probabilities are now given by, 



p(a'\<f>) = l/ntr[\a')(a'\pp ia (<f>)] 



lot' -a cxp[i<ft] | 

e -3+1 



n^/np + 1 



(58) 

From this the Fisher Information and the minimum vari- 
ance can be readily computed: 



np + 1 _ np + 1 
2iV|a| 2 ~~ 2Nn n 



(59) 



which again is sub-optimal for pure states (although only 
by a constant factor of 2) and approaches the optimal 
bound (|4"2")) at high temperatures. 



B. Canonical phase-measurement 

In the one-copy case we saw that the phase measure- 
ment, defined after Eq. ([5]), was optimal. Although some- 
what a theoretical exercise, it is interesting to study the 
accuracy of such measurement in the multicopy scenario. 
For pure squeezed states it is known that, although it is 
sub-optimal, the scaling of the variance with the mean 
photon number is the same as that for the optimal pro- 
tocol: Var[</>] ~ l/(Nn 2 ). Here we would like to study 
its performance in the presence of losses. The computa- 
tion for arbitrary squeezing is quite involved, but we can 
obtain analytical expressions in the two limiting cases 
of large and low initial squeezing. From Eq. (|29|) it is 
straightforward to get the expression of the probability, 



= £ (t 

n,m— 
n(n,m) 



n+m 



k=0 



where (f> = <fr — 8 and p^(k) are defined in (|5D)l . For 
large squeezing the sum above is dominated by large val- 
ues of n and m. Then the binomial distribution Pj^(k) 
is well approximated by a Gaussian distribution centered 
at 2nT and with variance 2nTR, and the sum over k can 
be easily performed as an integral (for that we can ex- 
tend the limits of integration to be ±oo). The remaining 
factorials can be approximated using the Stirling formula 



a.s 



,f{2n)\/n\ ~ 2™(mr)- 1 /4_ We obtain 



v/^A 2 



2tt 



E 2i(f>(n—m) \n+m 
e a 



(m + n)ir 



(61) 
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At this point it is convenient to define the variables s 
m + n and u = in — n and rewrite (1611) as 



p(4>) 



2jt 



s=0 J - 



du e 2s « e 



E('- ,5? *.)' 

s=0 



(62) 



where we have considered u as a continuous variable 
and have extended the integration range from (— s, s) to 
(— oo, oo). This is legitimated as it amounts to neglecting 
terms that fall off exponentially in u. The last sum can 
be trivially performed and yields 



p(4>) 



AovT 



X 2 



An 



(63) 



Now we can compute the Fisher information from 
Eq. (|40|) by noticing that for large squeezing (Ao — > 1) 
the angular integral is dominated by small values of cf> 
and, therefore, exp[2(T/i?)^ 2 ] ~ 1 + 2(T/R)4> 2 . Extend- 
ing to ±oo the range of the integral over 9, and recalling 
that in this limit Ao — 1 — l/(2n ), the Fisher information 
finally reads / ~ 2(T/R)n . Thus, the variance is 



Var-[0] 



1 — T 
2NTn~ Q 



(64) 



This value is again twice as large as the optimal variance 
in (|46[). and coincides with that of the heterodyne scheme. 

We next focus on the phase measurements in the op- 
posite limit: no < 1. Here the computation is much 
less involved. One just needs to keep terms up to or- 
der A 2 , in Eq. (|60|) , compute the Fisher information, and 
take its linearized expression to the same order Ag. Af- 
ter performing the trivial angular integration, we obtain 
/ ~ 4T 2 Aq ~ 4T 2 n (recall that in this limit A ~ y/m). 
Therefore the variance is 



Var" 



1 



ANT 2 n ' 



(65) 



Note that this value is twice the optimal variance in the 
same limit (|47|) . 

For coherent thermal states we proceed in a similar 
fashion. We first compute the phase-measurement out- 
come probabilities starting from the representation of 
displaced thermal states (fTTf for high-field amplitudes 
(a > 1), 



p(4>) 



l 

irnp 



da' 2 e ~e \S(Ce lv )\ 2 



(66) 



where we have implicitly defined £e iV = (a + o/)e l< ^, and 



where 



^{2-kC, 2 )^ 
2C 2 



ni 

oo 



due 



72 



yz a ■ yz 



(67) 



In the second equality we have used the Gaussian ap- 
proximation to the Poissonian photon distribution of a 
coherent state and extended the lower limit of integra- 
tion from zero to — oo. In the limit of large amplitudes 
£ f» a S> 1, and, hence, the distribution \S(£e llp )\ 2 is 
strongly peaked at <p w 0. Thus, the following approxi- 
mation holds 



(68) 



(ip w Csin^ = Im[(a + a')e^] 
= (a + a' x ) sin <j> + a' y cos <f>. 



With this we can and easily perform the integral in ([66 
since it becomes Gaussian. We obtain 



P$) 



2a 2 



■n{2n (i + 1) 



(69) 



It is now is straightforward to perform the 4> integral 
in the definition of the Fisher information, Eq. (|4"0"|) , and 
obtain the variance for the phasc-covariant measurement 



2n p + 1 _ 2n fj + 1 



(70) 



47V|a| 2 ANn a ' 
which agrees with the optimal bound. 

C. Homodyne measurements 

For pure states the standard quantum limit is known 
to be attained by homodyne measurements with a sim- 
ple two step adaptive protocol [l8[ . These measurements 
are very relevant from the practical point of view since 
they can be readily implemented in optical experiments 
[20l . [4l| . It seems natural to expect that homodyne mea- 
surements are also asymptotically optimal for thermal 
states, i.e., that they saturate the bound (|43[) . Let us 
anticipate that this is not the case. 

In the limit of strong local oscillator, a homodyne mea- 
surement is described by the set of projectors {|a;) e (a;|}, 
where are the eigenstates of the quadrature op- 

erator xg = (ae l6 + at e~ ,e ) / 'y/2. The probability of ob- 
taining outcome x upon measuring a shifted squeezed 
thermal state pp, r (<fi) is given by [42| . 



p g (x\<t>) = {x\U{4>-6)S{r)ppS{ryU{4>-ey\x) 

c 2 

■ ex P[-WZT77\ In]. ( 71 ) 



1 



^2tto- 2 ((9-0) 



2a 2 { 
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where 



a 2 {9) = (2np + l){e 



2r . 



e- 2r sin 2 



9). (72) 



The Fisher information can be easily computed from 
its definition in Eq. (|40|) . One obtains 



Ie(<t>) 



1 



sin[2(0- 0)](e 2r - e~ 2r ) 



(73) 



At this point it is readily apparent that the tempera- 
ture dependence has disappeared from the expression of 
the Fisher information. The maximum of the Ig(4>) is 
achieved for an homodyning angle 

6* = 0±arctan(e 2r ), (74) 

and the maximum Fisher information reads 



I* 



2sinh 2 (2r) = 8n r {n r + 1). 



(75) 



Upon using the Cramer-Rao bound we obtain that the 
maximum precision that can be achieved with homodyn- 
ing measurements is 



Var, 



horn r 

sq [ 



i 



8Nn r (n r 



1) 



(76) 



The above maximum of the Fisher information can only 
be attained if some previous information on the phase is 
available. It can be argued [l8[ that one can use a vanish- 
ingly small fraction of copies to perform a (perhaps) sub- 
optimal rough estimation of the phase and then use the 
remaining number of copies to nail down its precise value. 
The final estimation fidelity is dominated by this second 
stage, and therefore the above bound can be asymptot- 
ically attained. A similar analysis leading to the same 
results was recently carried out by Mitchell et qLp3|. 

Several comments are in order when comparing (|76|) 
with the optimal variance (|43|) . First, we observe that 
homodyning protocols are sub-optimal, they yield the 
optimal variance only for pure states. Second, we see 
that the variance (|76p is independent of the tempera- 
ture of the thermal state. We can understand this result 
as follows. For squeezed states, the effect of tempera- 
ture is the same as rescaling the phase-space. At the 
same time, homodyne measurements can be understood 
in phase-space as projections of the distribution on one 
quadrature. Therefore, the effect of temperature in the 
measurement statistics can be accounted for by a simple 
rescaling, thus the final precision is independent of the 
temperature. 

The variance ([76^1 for a squeezed state that emerges 
from a lossy channel of transmittance T can be calculated 
by using 



Var* 



l + 4n T(l-T) 
8NT 2 n (n Q + 1)' 



(77) 



where we recall that n = sinh r is the mean pho- 
ton number of the initial squeezed vacuum that is sent 



through the channel. Of course, for T = 1 one recovers 
the value of the variance for pure states [HI, For 
large squeezing, no 1, one has 



l-T 

2NTn. 



o 



while for no -C 1, 



Var; 



horn r 



l 

8NT 2 n ' 



(78) 



(79) 



So, in the presence of losses and for low squeezing, adap- 
tive homodyning provides a variance a factor of 2 smaller 
than the variance of heterodyning, and attains the opti- 
mal bound (|46p . while for high squeezing homodyning 
provides half the precision of the optimal strategy but 
still outperforms heterodyning by a factor which depends 
on the amount of losses. 

For coherent thermal states, the same scaling argument 
we used for squeezed states does not hold because the dis- 
placement is not affected by an increase in temperature, 
while it is obviously affected by phase-space rescaling. It 
is however a straightforward exercise to check that the 
Fisher information for displaced thermal states is 



H4>) 



4|a| 2 sin 2 ( 



2np + 1 



(80) 



Therefore, the optimal variance (|42|) for displaced ther- 
mal states can be attained with a one-step adaptive ho- 
modyne measurement by taking 6* = <p ± 7r/2, even for 
states with non-zero temperature, 



Var c h o ™ = 



1 2-np + l 
AN n n 



(81) 



We notice here that the Fisher information, or the min- 
imum variance for that matter, are the same for a 
coherent state \a) than for a displaced thermal state 
D(a')ppD(a') f with \a'\ 2 = \a\ 2 (2np + 1), which can be 
understood in terms of the scaling argument used above. 

It still remains an open problem to find optimal and 
practical strategies that attain the bound for general 
squeezed thermal states. We note that the performance 
of adaptive homodyne measurements is as good as, or 
better, than that of the single-shot canonical phase mea- 
surement for both displaced and squeezed thermal states. 



VI. FREQUENCY ESTIMATION 

Sub shot-noise parameter estimation is one of the most 
important applications brought by the field of quantum 
information. Since the work of Huelga et al. [H it is well- 
known that the use of exotic quantum states provides 
a dramatic improvement on the estimation precision of 
several metrologically relevant parameters. It was imme- 
diately recognized that states that provide this enhance- 
ment can be specially fragile in the presence of decoher- 
ence or losses and that it is therefore mandatory to asses 
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whether the quantum advantage persists in noisy envi- 
ronments @, [3] • Recently, there have been several in- 
depth studies on this issue in discrete systems, specially 
oriented to the so-called maximally path-entanglement 
state (NOON) state (\N,0) + |0,^»g|liSSl- 

Here we want to approach the problem for CV systems. 
In particular we will use the previous minimum values for 
the phase estimation variance to study a model for fre- 
quency estimation in the presence of losses (which is the 
dominant source of errors in optical experiments). We 
will study a scenario where the phase is imprinted at a 
given rate <p = ut on a Gaussian state, in a time interval 
t during which the system suffers losses also at a con- 
stant rate r\. Long times t will increase the accumulated 
phase and hence improve the sensitivity of the measure- 
ment. However, the unavoidable losses will clearly limit 
the duration of the experiment. Our goal here is to find 
the optimal time t* to estimate the frequency u with the 
highest precision. 

As stated above, a BS can be used to model a lossy 
channel. In our scenario the input state, a pure squeezed 
or displaced state, is sent through a BS with a time de- 
pendent transmittancc T, T = e~ r,t . So, we can use the 
above results, (|4"2"| and to optimize the time by tak- 
ing into account that now the variance for a fixed number 
of copies TV is given by 



Var[w] = [Nt 2 I\ 



(82) 



The case of coherent states can be solved analytically. 
When a coherent state |a) is sent to a BS, the transmit- 
ted state is nothing but \T x / 2 a). So, we only have to 
optimize over time the following expression: 



Var coh [oj] 



1 



AN\a\ 2 e-^t 2 



(83) 



where we have replaced \a\ 2 by \a\ 2 e~ vt and taken the 
limit (3 — > 00 in Eq. (|4"2"|) . We find that the optimal time 
is t* = 2/r), which gives the minimum variance 



Var coh [w] 



2 2 
e 77 

16N\a\ 2 ' 



(84) 



For squeezed input states we recall that the squeezing 
parameter and the temperature of the state at a given 
time t depend on T(t) = e -11 and on the initial squeez- 
ing ro through Eq. (|28[) . Inserting this time dependence 
in (|4"5")) we can optimize Var sq [w] = Var sq [0]/t 2 . Fig- 
ure [3] shows the optimal value of rjt* for different val- 
ues of ro together with the corresponding minimum vari- 
ance. The optimal time is now a function of the initial 
squeezing and ranges from t* = l/rj for low squeezing 
to t* = (2 + W(— 2e~ 2 ))/r] f» 1.59/7? for high squeezing, 
where W(x) is Lambert's TF-function [Zj|. At this high 
squeezing regime the variance scales as 0(n,Q 1 ). Hence, 
the characteristic Hciscnbcrg scaling of squeezed pure 
states turns into the standard quantum limit scaling in 
the presence of losses. 



0.01 




FIG. 3: (Color online) Rescaled variance iV»7 2 Var[a;] (solid 
line) and optimal time r\t* (dashed line) for the case of input 
squeezed vacuum state with squeezing parameter ro. 
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FIG. 4: (Color online) Rescaled optimal variance u>\ 
for coherent (solid) and squeezed (dashed) states as a function 
of the mean photon number no of the initial state. The dotted 
line is the squeezed state variance for an adaptive homodyne 
measurement, while the long-dashed line is the variance for 
heterodyne measurements. Inset: As above but for a range of 
small values of no. 



In Figure [5] we compare the performance of an initial 
pure coherent state with that of a pure squeezed state of 
the same energy (mean photon number). In the limit of 
low photon numbers (sec inset) the squeezed state per- 
forms worse than the corresponding coherent state of the 
same energy. At higher energies the squeezed state gives 
a slightly higher precision. 

A similar behavior is found if one restricts to the adap- 
tive homodyning strategy mentioned above. In that case 
the variance in presence of losses for a given time t can be 
computed from (|77|) . The optimal times arc shorter than 
those required for the optimal POVM, but they have the 
same asymptotic values r\t* G [1, (2 + W(-2e~ 2 ))]. Fig- 
ureUshows the corresponding minimum variances for the 
homodyning strategy. We note that for small squeezing, 
homodyning is nearly optimal, while for large squeezing 
the differences become important (close to a factor of 2), 
though still give rise to the same type of scaling. 

Finally, using (|5"B"|) we can obtain the results for the 
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heterodyne strategy. In this case the optimal times 
range from rjt* = 1 for very low squeezing to rjt* = 
2 + W{— e~ 2 ) fa 1.84. As shown in FigurclH heterodyn- 
ing always performs worse than homodyning. The dif- 
ferences are especially large in the low squeezing regime, 
where homodyning is nearly optimal. 

VII. CONCLUSIONS 

We have derived optimal bounds to the phase estima- 
tion fidelity for pure and mixed Gaussian states of light. 
This enables us to study how temperature affects the pre- 
cision of the estimation. We have focused on displaced 
thermal and squeezed thermal states. A priori, one would 
expect that the fidelity degrades with temperature, as the 
fluctuations also increase. This is precisely what we find 
for coherent thermal states. However, squeezed thermal 
states have the opposite behavior: a squeezed state with 
high temperature can provide twice the precision of a 
pure squeezed state with the same squeezing parameter. 
The result seems at first sight paradoxical, but it can be 
understood by noticing that an initial thermal state gives 
rise to a non-linear increase in the mean photon number, 
which counterbalances the thermal noise. 

In the many-copy scenario we have studied three differ- 
ent single-shot measurements: heterodyning, canonical 
single-copy phase measurement, and adaptive homodyn- 
ing. Adaptive homodyne measurements attain the op- 
timal bounds for coherent thermal states and squeezed 
vacuum states but are sub-optimal for thermal squeezed 
states since they provide a precision that is independent 
of temperature. Heterodyning provides sub-optimal pre- 
cision for pure states; however we find the surprising 
result that it does saturate the optimal bounds in the 



limit of very mixed states, a regime where adaptive ho- 
modyning ceases to be optimal. Finally, what is taken 
to be the canonical phase measurement, which is opti- 
mal in the single-copy scenario but extremely hard to 
implement, turns out to perform worse than the practi- 
cal heterodyne or adaptive homodyne strategies in some 
cases when many copies are available. 

We have given a simple model to study the effect of 
losses in a situation where the phase imprinted on a gaus- 
sian state grows linearly with time, <f> = ut, while the 
system suffers losses also at constant rate rj. We give the 
optimal duration of the experiment which finds the com- 
promise between losses and accumulated phase. Using 
the derived many-copy results and optimizing over this 
time interval, we show that high-energy squeezed states 
give only a minor advantage over coherent states with the 
same energy, and are in fact outperformed by them in the 
low energy regime. While a similar analysis could be car- 
ried out for other Gaussian noisy channels (the complete 
Bures metric in the Gaussian state space, from which 
one obtains the optimal bound to the variance, can be 
found in [48[), it remains an open problem to address op- 
timal estimation for non-Gaussian noisy channels, such 
as phase diffusion. 
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